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ABSTRACT 

Absorbers seen in the spectrum of background quasars are a unique tool to select 
Hl-rich galaxies at all redshifts. In turns, these allow to determine the cosmological 
evolution of the H i gas, 57hi+HoII, a possible indicator of gas consumption as star 
formation proceeds. The Damped Lyman-a systems (DLAs with A^hi^ lO^" '^ cm~^), 
in particular, are believed to contain a large fraction of the H i gas but there are 
also indications that lower column density systems, named "sub-Damped Lyman-a" 
systems play a role at high-redshift. Here we present the discovery of high-redshift 
sub-DLAs based on 17 z > 4 quasar spectra observed with the Ultraviolet- Visual 
Echelle Spectrograph (UVES) on VLT. This sample is composed of 21 new sub-DLAs 
which, together with another 10 systems from previous ESO archive studies, make up 
a homogeneous sample. The redshift evolution of the number density of several classes 
of absorbers is derived and shows that all systems seem to be evolving in the redshift 
range from z = 5 to z ~ 3. This is further used to estimate the redshift evolution of 
the characteristic radius of these classes of absorbers assuming a Holmberg relation 
and one unique underlying parent population. DLAs are found to have i?* ^ 20ft-]^oo 
kpc, while sub-DLAs have i?* ~ 40/ij~QQ kpc. The redshift evolution of the column 
density distribution, f(N,z), down to Nm= 10^^ cm~^ is also presented. A departure 
from a power law due to a flattening of f(N,z) in the sub-DLA regime is present in the 
data. f(N,z) is further used to determine the H i gas mass contained in sub-DLAs at 
z>2. The complete sample shows that sub-DLAs are important at all redshifts from 
z = 5 to z = 2. Finally, the possibility that sub-DLAs are less affected by the effects 
of dust obscuration than classical DLAs are discussed. 
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1 INTRODUCTION 

Tracing the rate at which stars form over cosmological scales 
still remains a challenging observational task. An indirect 
way to probe the assembly of galaxies is to probe the rate 
at which they convert their gas into stars. The neutral H I 
mass in particular can be estimated from observations of ab- 
sorbers seen in the spectrum of background quasars. Unlike 

* Based on observations collected during programme ESO 69. A- 
0613, ESO 71.A-0114 and ESO 73.A-0653 at the European South- 
ern Observatory with UVES on the 8.2 m KUEYEN telescope 
operated at the Paranal Observatory, Chile, 
t e-mail: cperoux@eso.org 



other high-redshift galaxies (such as Lyman Break Galaxies, 
Steidel et al. 2003), these objects are selected solely on their 
H I cross-sections regardless of their intrinsic luminosities or 
star formation rates. The quasar absorption systems are di- 
vided into several classes according to the number of atoms 
along the observed line of sight: the Lya forest have H I col- 
umn densities ranging from ~ lO^'^ to 1.6x 10^^ atoms cm"'^, 
the Lyman-limit systems (LLSs) with N(HI) > L6 x lO" 
atoms cm~^ and the damped Lya systems (DLAs) with 
N(HI) > 2 X 10^" atoms cm~^. These latter systems are be- 
lieved to be the major contributors to the neutral gas in the 
Universe at high redshifts. They thus can be used to mea- 
sure the redshift evolution of nni+Heii, the total amount of 
neutral gas expressed as a fraction of today's critical density 
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(Lanzetta et al. 1991; Wolfe et al. 1995; Storrie-Lombardi et 
al. 1996a; Storrie-Lombardi & Wolfe 2000). We have recently 
suggested (Peroux, et al. 2003a) that at z > 3.5, some frac- 
tion of the H I lies in systems below the traditional DLA 
definition, in "sub-Dampod Lyman-Q Systems (sub-DLAs)" 
with 19 < log N(HI) < 20.3 cm~'^. The name arises from the 
fact that this column density range lies on the linear part 
of the curve of growth which links the equivalent width of 
an absorber to the number of atoms along the line of sight. 
The present paper reviews these predictions based on direct 
observations of a sample of sub-DLAs. 



In addition, quasar absorbers arc a direct probe of 
element abundances over > 90% of the ago of the Uni- 
verse. The cosmological evolution of the HI column den- 
sity weighted metallicity (e.g., Kulkarni & Fall 2002) shows 
surprising results: contrary to virtually all chemical models 
(e.g., Malanoy & Chaboyor 1996; Pei, Fall & Hausor 1999), 
the most recent observations indicate mild evolution with 
redshift (Khare, et al. 2004; Kulkarni, et al. 2005). There- 
fore, the results to date might give a biased or incomplete 
view of the global galactic chemical evolution. The ultimate 
goal of this project is to undertake detailed abundance and 
dynamical studies of an appropriate sample of sub-DLAs 
in order to estimate the evolution of their mean abundance 
with redshift. 



In a first step towards this goal, we took advantage of 
the ESO VLT archive to build a sample of sub-DLAs by 
analysing UVES archival echelle QSO spectra. This repre- 
sents a sample of 35 QSOs, 22 of which were unbiased for 
our study. This study led to the discovery of 12 sub-DLAs 
(Dessauges-Zavadsky et al. 2003). Their chemical abun- 
dances were derived using Voigt profile fitting and photoion- 
ization models from the CLOUDY software package in order 
to determine the ionization correction. We find that the cor- 
rection is negligible in systems with N{HI) > 3.2 x 10^^ and 
lower than 0.3 dex for most element in systems with 10^^ < 
N(HI) < 3.2 X lO^'' atoms cm~^. These systems were used to 
observationally determine the shape of the column density 
distribution, f{N,z), down to N(HI) = 10^® cm"^ (Peroux 
et al. 2003b), although the lack of high-redshift systems pre- 
vented from measuring the f{N, z) redshift evolution. The 
abundances observed in this sample of sub-DLAs were fur- 
ther used to determine the global metallicity of H i gas in 
both DLAs and sub-DLAs. 



Here we present a new sample of 17 z>4 quasar lines of 
sight observed at high-resolution with the Ultraviolet- Visual 
Echelle Spectrograph (UVES) on VLT. We use these newly 
acquired data to search for and study the statistical proper- 
ties of high-redshift sub-DLAs. In the second section, the ob- 
servations and data reduction process are detailed together 
with a description of each object. The third section presents 
the statistical characteristics of the survey and the newly 
built sub-DLA sample. Finally, our current state of knowl- 
edge of this class of absorbers together with implications for 
the cosmological evolution of H I gas mass are discussed in 
section 4. 



2 THE Z >2 SUB-DLA SAMPLE 

2.1 Observations and Data Reduction 

In order to complement the previous study and analyse sub- 
DLAs at z > 3, we started in September 2002 to build a 
UVES sample of high redshift QSOs never observed before 
at high resolution. This snapshot is carefully designed to 
determine the redshift evolution of the statistical properties 
of sub-DLAs (rmmbcr density and column density distribu- 
tion). The second step of the observing programme under- 
took a year later, in September 2003, aims at concentrating 
on the 7 most promising targets to enable detailed metal- 
licity and dynamical studies at high-redshifts. An extra set 
of data was obtained on August 29"" 2004, September 2""^ 
and 4"^ 2004 in service mode in the framework of a paral- 
lel observing programme (ESO 73.A-0653; PI: N. Bouche). 
Table 1 gives the journal of the observations. 

The data reduction of the observations from period 
69 and 71 is done using the version of the UVES pipeline 
within the MIDAS environment available at the time (ver- 
sion: uves/2.0.0 flmidas/l.O.O). The newest data (period 73) 
are reduced with the most recent version of the pipeline to 
accommodate for new format of the raw fits file (version: 
uves/2.1.0 flmidEis/ 1.1.0). The raw wavelength frames are 
inspected closely to optimise the number of order to be ex- 
tracted for each CCD. Master bias and fiat images axe con- 
structed using calibration frames taken the closest in time 
to the science frames. Lamps arc taken on various occasions 
through the observing nights to be able to make accurate 
wavelength calibrations. In most cases, the science frames 
and associated stars are extracted with the "optimal" op- 
tion. Nevertheless, on a few frames from the second observ- 
ing run (September 2003), two of the settings (540 and 800) 
show a periodic feature of the order of a few Angstroms. 
These features are not present when these few frames are 
extracted using the "average" option followed by a cosmic 
filtering. 

For the 7 objects with more than one setting, the re- 
sulting spectra are combined by using the signal-to-noise of 
the spectrum as a weight. The final spectra have resolution 
of 7 km s~^. In some case, we notice that the part of the 
spectra fully absorbed do not actually reach the zero flux 
level. This problem is known to occur when there axe not 
enough photons in the sky area. Indeed, the optimal extrac- 
tion is done to fit a GauBian profile and the baseline of the 
Gaufiian is the sky level. If the seeing during the observa- 
tions is poor and the exposure short, it is possible that at 
blue wavelengths, only a few photons fall on the sky pixel. 
This leads to a high photon noise and might produce the 
under-subtraction of the sky level. The unsaturated absorp- 
tion lines are unaffected by this zero-level problem. In the 
case of saturated lines, we correct the spectrum by subtract- 
ing a few percent from the continuum level. The spectra are 
then normalised using a spline function to join the part of 
the continuum free from absorption lines. 

2.2 Absorbers Identification and Column Density 

Measurements 

The methodology to search for sub-DLAs in the spectra de- 
scribed above very much follows the one from our previ- 
ous studies (Dessauges-Zavadsky et al. 2003). We first use 
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Table 1. Journal of observations of our observing programmes to find and analyse new high-redshift sub-DLAs in 17 2 > 4 quasars. 



Quasar 


Alternative name 




R mag 


Obs Date 


UVES settings 


Exp. Time (sec) 


Ref 


hSK JUUUo— DzUo 




4.455 


19.25 


z( ocpt ZUUz 


580 


7200 


1 


Fbo JUllo+UozU 




4.230 


18.66 


2o bcpt 2[)02 


560 


3600 


6 










Ul oept zUUo 


oyO+580 


5400 












Ul oept zUUo 


580 


5400 












A1 CJj^r^-t- OAAQ 

Ul oept zUUo 


860 


04UU X Z 




Foo JU12i+Uo47 




A 1 0*7 


18.33 


OO CJ^T^+ OAAO 

2o oept zUU2 


540 


3600 


5 










Uz bept 2UUo 


47U+oUU 


JdUUxz 




oUbb JU124+UU44 




3.840 


17.75 


OV C^i*-.+ OAAO 

2/ bept ZUUz 


520 


3600 


A 

4 










AO C^i^+ OAAQ 

Uz bcpt zUUo 


520 


3600 












OA /A/I a .1^1- /C^*-.+ OAA/1 

zy/U4 Aug/ bept ZUU4 


580 


Q AAA \f A 

oUUUX4 




Fbb JUloo+U4UU 




4.154 


18.32 


Z7 bcpt zUUZ 


520 


3600 




"DTDT 7A1 O'? Ann Aved 
ortl JUlo/— 4zz4 


ortl rSUioO— 4zoy 


3.970 


18.77 


OT OAAO 

z7 oept zUUz 


520 


3600 


3 










AQ CJj^r^-i- OAAQ 

Uo oept zUUo 


COA 

ozU 


5000 












AQ O^r^T^-i- OAAQ 

Uo oept zUUo 


A TA 1 OAA 

47UH-oUU 


QfiAA 1 A CAA 

ooUU-f-4oUU 




occ TAOAn 1 1 '7 
Foo JUzUy+UolY 




A t TA 
4.1/4 


17. YD 


OO C!^T^+ OAAO 

2o oept zUUz 


QOA 1 CQA 

oyU+ooU 


7200 


1 T 


rSHi JUz44— Uio4 


rSrti JDUz4i — U14D 


4.053 


18.18 


OT O^i^-f- OAAO 

Z ( bcpt ZUUZ 


520 


3600 


2 


"DTD TAQI 1 1 '700 




4.039 


17.73 


OO Gr^nA- OAAO 

Zo bcpt ZUUZ 


540 


3600 


1 






A Q^iQ 
4.000 




OC G«i-,-(- OAAO 

Zo bcpt ZUUZ 


OOU 


QfiAA 
v>DUU 


± 










AO Gr^nA- OAAQ 

UZ bcpt ZUUo 


580 


obUUXZ 












AO Gr^nA- OAAQ 

UZ bcpt ZUUvS 


860 


5066 




■O T> Tr\/11A CTlii; 

BH JU4iy— o71d 




4.461 


17.78 


OO O 4- OAAO 

2S oept zUUZ 


560 


1846 


1 




XjXV-Li v_-' w JjZUIO — "iUZo 




ICS. ou 


Ut> oepx ZUUo 


04U 


OOUU 


1 


PSS J2154+0335 




4.363 


19.05 


28 Sept 2002 


560 


3600 


1 


BR J2215-16ir^'^ 


BR B2212-1626 


3.990 


18.10 


28 Sept 2002 


540 


3600 


2 










03 Sept 2003 


540 


3600+1215 












03 Sept 2003 


470+800 


3600x2 




BR J2216-6714''^'^ 




4.469 


18.57 


28 Sept 2002 


560 


3600 


1 










01 Sept 2003 


560 


3600x2 












01 Sept 2003 


580 


5400 












02 Sept 2003 


860 


5400x2 




BR J2239-0552 


BR B2237-0607 


4.558 


18.30 


27 Sept 2002 


580 


3600 


2 


BR J2349-3712 


BRLCO B2346-3729 


4.208 


18.70 


03 Sept 2003 


540 


3600 


1 



BAL. affected by Broad Absorption Lines features. 

spectrum with complete wavelength coverage including to the red of the Ly-a emission line. 
^<''^*: attempt to get a spectrum with complete wavelength coverage impaired by bad weather conditions 
(strong northern wind). 

References - (1) Pcroux ct al. 2001; (2) Storric-Lombardi ct al. 1996a; (3) Storrie-Lombardi et al. 
2001; (4) Schneider ct al. 2002; (5) Stern et al. 2000; (6) from the Palomar Sky Survey, see 
http://www.astro.caltech.edu/~george/z4.qsos; (7) Prochaska et al. 2003. 



an automated detection algorithm supplemented further by 
eye-baJling searches independently undertook by three of us 
(C. P., M. D.-Z. and T.-S. K.). Although the forest of the 
quasars at these redshifts is considerably absorbed by fore- 
ground structures, the presence of damping wings down to 
the Ly2 and occasionally to higher lines of the series provides 
an unambiguous signature of damped absorbers. We also re- 
port the detection of metal lines at the associated redshifts, 
although this is not used as a criterion for sub-DLA selec- 
tion. Table 2 summarises these information together with 
the redshifts, H i column densities and the references of the 
7 DLAs and 21 sub-DLAs found in our new high-redshift 
quasar sample. The latter systems span a range of H I col- 
umn densities from lO^^ "^^ to lO^''-^^ cm"^ with 2abs=2.976 
to 4.145. 

The H I column density axe determined by fitting a 
Voigt profile to the absorption line. The fits were per- 



formed using the minimisation routine f itlyman in MIDAS 
(Fontana & Ballester 1995). The Doppler parameter 6- value 
was usually fixed at 20 km or left as a free parame- 
ter since in that high column density regime, the sub-DLA 
systems Nbi are indepeii<leiit of the 6-value. The fit is per- 
formed using the higher members of the Lyman series where 
these are available. The typical resulting error bar in log iVni 
measurements is 0.10 and never exceeds 0.15, but does not 
include error in the coritirmum placement which we expect 
should not exceed 10%. 



2.3 Notes on Individual Objects 

In this section, we provide details on each of the identified 
absorbers. Whenever the spectral coverage is available, we 
search for metals associated with the Ly lines among the 20 
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Table 2. High-redshift quasar absorber sample composed of 21 sub-DLAs and 7 DLAs. 



Quasar 




■^abs 


log A^Hi 


Ly scries 


metals 


Note 


BR J0006-6208 


4.455 


3.202 


20.80±0 


10 


Lyl 


no inetals over available coverage 


DLA 






4.145 


19.37±0 


15 


Ly3 


no metals over available coverage 




PSS J0118+0320'"<='' 


4.230 


4.128 


20.02±0 


15 


Ly5 


Fell, Sill, OI, CII, SilV, CIV 




PSS J0121+0347''«^* 


4.127 


2.976 


19.53±0 


10 


Lyl 


OI, CIV, Fell, AlII 




SDSS J0124+0044''«<^ 


3.840 


2.988 


19.18±0 


10 


Lyl 


Sill, CIV, Fell, AlII 








3.078 


20.21±0 


10 


Lyl 


CII, Sill, OI, SilV, CIV, Foil 




PSS J0133+0400 


4.154 


3.139 


19.01±0 


10 


Lyl 


Nill, CII 








3.692 


20.68±0 


15 


Ly2 


no metals over available coverage 


DLA 






3.773 


20.42±0 


10 


Ly2 


no metals over available coverage 


DLA 






3.995 


19.94±0 


15 


Ly4 


no metals over available coverage 


blended with the following system 






3.999 


19.16±0 


15 


Ly4 


no metals over available coverage 


blended with the previous system 






4.021 


19.09±0 


15 


Ly4 


Fell 




BRI J0137-4224''«<^ 


3.970 


3.101 


19.81±0 


10 


Ly4 


Sill, CIV, Fell, AlII 








3.665 


19.11±0 


10 


Lyll 


Sill, OI, CII, SilV 




PSS J0209+0517 


4.174 


3.666 


20.47±0 


10 


Ly2 


SilV 


DLA 






3.707 


19.24±0 


10 


Ly2 


01, SilV 








3.863 


20.43±0 


15 


Ly3 


Sill, OI, CII 


DLA 


BRI J0244-0134 


4.053 














BR J0311-1722 


4.039 


3.734 


1 o /I c_i_n 


1 fl 


Ly7 


Sill, OI, CII 




BR J0334-1612'''^'' 


4.363 


3.557 


21.12±0 


15 


Lyl 


Sill, CIV, Fell, AlII, Znll 


DLA 


BR J0419-5716 


4.461 


3.063 


19.17±0 


10 


Lyl 


OI 




BR J2017-4019^'*^ 


4.131 














PSS J2154+0335 


4.363 


3.177 


19.23±0 


15 


Lyl 


no metals over available coverage 




BR 12215-1611''='* 


3.990 


3.656 


ig.oiito 


15 


Ly3 


Nill, CII, Znll 


blended with the following system 






3.662 


20.05±0 


15 


Ly3 


OI, CII, SilV, Sill, CIV, Fell, AlII 


blended with the previous system 


BR J2216-6714'''^'' 


4.469 


3.368 


19.80±0 


10 


Lyl 


CII, CIV, AlII 




BR J2239-0552 


4.558 


4.079 


20.55±0 


10 


Ly3 


Sill, OI, CII 


DLA 


BR J2349-3712 


4.208 


3.581 


19.12±0 


10 


Ly2 


no metals over available coverage 








3.690 


19.79±0 


15 


Ly2 


no metals over available coverage 


blended with the following system 






3.696 


19.78±0 


10 


Ly2 


no metals over available coverage 


blended with the previous system 



BAL. affected by Broad Absorption Lines features. 

spectrum with complete wavelength coverage including to the red of the Ly-a emission line. 
^<"i*: attempt to get a spectrum with complete wavelength coverage impaired by bad weather conditions 
(strong northern wind). 



transitions most frequently detected in high column density 
quasar absorbers. 

(i) BR J0006-6208 («em=4.455): Peroux et al. (2001) 
have reported four candidate absorbers at «abs=2.97, 3.20, 
3.78 and 4.14 in the medium resolution spectrum of this 
qucisar. The blue end of the spectrum presented here is 
of fairly poor quality and therefore the minimum redsliift 
above which quasar absorbers are searched for is set to 
Zmin=3.083, so the lowest DLA cannot be checked upon. 
We confirm the presence of the 2:abs=3.202 system and mea- 
sure log A'^Hi=20.80±0.10 in agreement with medium reso- 
lution estimates. The 3rd system falls in the UVES set- 
ting gap. The last absorber is confirmed as a sub-DLA with 
logiVHi=19.37±0.15. No associated metal lines are detected 
over the limited wavelength coverage of our spectrum. Fig- 
ure 1 shows the best Voigt profile fitting solutions for these 
two systems. 

(ii) PSS J0118-F0320 (aem=4.230): This quasar was dis- 
covered as part of the Palomar Sky Survey and to our knowl- 
edge, there are no medium resolution spectrum published. 



We discover an unambiguous sub-DLA at 2:abs=4.128 and 
measure log A'hi=20.02±0.15 down to the Lyman 5. Several 
associated metals are also observed in the red part of the 
spectrum. Figure 2 shows our best fit of the H i lines. 

(iii) PSS J012H-0347 (2em=4.127): No absorbers have 
been previously reported in this quasar also from the Palo- 
mar Sky Survey. We report here for the first time a sub-DLA 
with log AfHi = 19.53±0.10 at Zabs=2.976. Metal lines at that 
redshift are also observed in the red part of the spectrum. 
Figure 3 shows our best fit of the H I line. 

(iv) SDSS J0124-f0044 (zem=3.840): One absorber have 
been previously reported in this quasar from the Sloan 
Digital Sky Survey (Bouchc & Lowentlial, 2004). We re- 
port here for the first time the column density of the 
two sub-DLAs with logiVHi=19.18±0.10 at Zabs=2.988 and 
logiVHi=20.21±0.10 at 2:abs=3.078. Metal lines at these red- 
shifts are also observed in the red part of the spectrum. Fig- 
ure 4 shows our best fit of the H I lines. 
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BR J0006-6208 
DLA z=3.202 Lyi 




-2000 O 2000 

Relative Velocity (km/s) 

■ ■ ' BR ■j6o66-6268' 
Sub z=4.145 Ly3 



Cir BR J0006-6208 ^ 



Sub z=4.145 Ly2 




-500 500 

Relative Velocity (km/s) 



Figure 1. The 

two absorbers detected towards BR J0006— 6208 at z^]^s=S. 202 
(logAfHi=20.80±0.10) and 2abs=4.145 (log AfHi=19.37±0.15). 



(v) PSS J0133+0400 (zom=4.154): Peroux et al. (2001) 
report two DLAs in this quasar which were further con- 
firmed by ESI/Keck II observations of Prochaska et al. 
(2003). Here, we measure log 7Vhi=20.68±0.15 at 2abs=3.692 
and logiVHi=20.42±0.10 at 2abs=3.773 (Prochaska et al. 
measure log iVHi=20.70lE!:i° and logNm^iO.b^to.ll respec- 
tively). We also note that these two systems are close by 
each other: there are separated by ~ 5170 km s~^with 
an additional log A^hi< 19.0 in between them (see second 
panel of Figure 5). This group constitutes almost a mul- 
tiple DLA (Lopez et al. 2001, Ellison & Lopez 2001, 
Lopez & Ellison 2003). In addition, we report the dis- 
covery of a further 4 sub-DLAs along the same line 
of sight. We measure log A^hi=19.01±0.10 at 2;abs=3.139, 
logArHi=19.94±0.15 at 2abs=3.995, log AfHi=19.16±0.15 at 
2abs=3.999 and log AfHi = 19.09±0.15 at Zabs=4.021. Again, 
two of these are very close by 2:abs=3.995 and 3.999 corre- 
sponding to 250 km s~^. There are fitted together using Ly2 
and Ly4 since those provide better constrains than Lyi. The 
Lyman series down to Ly4 is also available for the last 3 sub- 
DLAs. The various Lyman series absorption lines along this 



. 1 1 1 1 1 1 


■ PSS J'01 18+0320 
Sub z=4.128 Ly5 


1 1 1 1 1 1 1 . 


. 1 1 I 1 1 1 

- i J7 i 

TT^TWf ' iyiii|iiii|'i'i 


■ PSS J01 18+0320 
Sub z=4.128 Ly4 




■ PSS J'01 18+0320 
Sub z=4.128 Ly3 










. 1 1 1 1 1 


■pSS JOl 18+0320' 
Sub z=4.128 Ly2 


' 1 ' ' ' 1 - 


.1 1 1 1 1 


■ PSS 30118+0320 
Sub z-4.128 Lyi 


1 1 1 1 1 1. 


, , , 1 




' 



Relative Velocity (km/s) 

Figure 2. The absorber detected towards PSS J0118-t-0320 st 
2abs=4.128 (logArHi=20.02±0.15). 



PSS JOl 21 +0347 
Sub z=2.976 Lyi 




-2000 -1000 1000 2000 

Relative Velocity (km/s) 

Figure 3. The absorber detected towards PSS J0121+0347 at 
2abs=2.976 (logAfHi=19.53±0.10). 



rich line of sight are displayed in Figure 5. Given the limited 
wavelength coverage of our spectrum, only metals associated 
with the lowest and the highest system are detected. 

(vi) BRl J0137-4224 (2cm=3.970): Storrie-Lombardi 
et al. (2001) do not find any DLA in the spectrum of 
this quasar. We confirm that there are none of the high- 
est H I column densities, but we find two sub-DLAs. The 
first system is at Zabs=3.101 and has a H I column density 
log A'hi=19.81±0.10, while the second is at 2abs=3.665 and 
has log AfHi=19.11±0.10. Both these systems have metals as- 
sociated with them. The lyman series down to Lyll is also 



© 0000 RAS, MNRAS 000, 1-18 



6 Celine Peroux et al. 



"5 



PSS J01 33+0400 
Sub z=3.139 Lyi 




-2000 -1000 



1000 200 



Relative Velocity (km/s) 



X 
3 

L. 
« 



PSS J01 33+0400 
DLA z»3.692 Ly2 




I I I I I I I I I I I I I I I I I I 

PSS J01 33+0400 
DLA z=3.692 Lyi 



PSS J01 33+0400 
DLA z=3.773 Ly2 





PSS J01 33+0400 
DLA z=3.773 Lyi 



-2000 -1000 1000 2000 

IMotiw VMod^ (km/s) 



-2000 -1000 1000 2000 

RdaUM VMod^ (km/*) 



t i^^lidiy+'o^db j t te^jbii3;o:H)6 

Sub z-3.999 Ly4 ' ■ Subz-4.021Ly4 

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 



. 



Sub z-4.021 Ly3 



^ |is^WiW+'o4ob « 

><!■ S,ub z-3.999 Ly3 H x 

, . . „ . ,^ _ , 

.2 Aok'i/t^WxMAA .2 » ^ /aAi!»W d 



Sub z-3.999 Ly2 



. 



^ rts'j6ii3+d4()o 

Sub z-3.999 Lyi 




PSS J0133+040CI ' 
Sub z-4.021 Ly2 



p^sJoii3+046o 

Sub z-4.021 Lyi 




-2000 -1000 1000 2000 

RekiOvaVWoGilyOim/i) 



-1000 1000 

RakiOvaVWoGi^Oiin/i) 



Figure 5. The rich hne of sight towards PSS J0133+0400 is composed of the systems with following redshifts: ^^abs — 3. 139 
(logArHi=19.01±0.10), 2abs=3.692 (log AfHi=20.68±0.15), 2abs=3.773 (log ArHi=20.42±0.10), 2abs=3.995 (log AfHi=19-94±0.15) together 
with 2abs=3.999 (log ArHi=19.16±0.15) and finally 2abs=4.021 (log AfHi=19.09±0.15). © 0000 RAS, MNRAS 000, 1-18 
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Figure 4. The absorbers detected towards SDSS J0124+0044 
at Zabs=2.988 (log AfHi=19.18±0.10) and Zabs=3.078 
(logAfHi=20.21±0.10). 



visible for the second system. Figure 6 shows the best Voigt 
profile fitting solutions for these two systems. 

(vii) PSS J0209+0517 (2om=4.174): Peroux et al. (2001) 
report two DLAs in this quasar which were further confirmed 
by ESI/Keck 11 observations of Prochaska et al. (2003). 
Here, we measure log A''hi=20.47±0.10 at 2:abs=3.666 and 
log A'hi=20.45±0.15 at Zabs=3.863 (Prochaska et al. mea- 
sure logAfHi=20.55±0.10 and log iVHi=20.55±0.10 respec- 
tively). We also find a sub-DLA at 2:abs=3.707. We measure 
log A'hi = 19.24±0.10 for that system. All these absorbers 
have metals associated with them. Figure 7 shows our best 
fit of these H I lines. 

(viii) BR J0244-0134 (zcm=4.053): This quasar was ob- 
served by Storrie-Lombardi et al. (1996b) who did not re- 
port any DLA. The high-resolution spectrum that we have 
acquired shows that no sub-DLA is found either along this 
line of sight. 

(be) BR J0311-1722 (2om=4.039): Peroux et al. (2001) 
report an absorber with column density below the classi- 
cal definition along this line of sight from medium reso- 
lution spectroscopy. Here, we confirm that the system has 
logiVHi = 19.48±0.10 at Zabs=3.734 observable down to Ly7 
as shown in Figure 8. Most of the constraints in this system 
comes from Ly3 and Ly5. The lower members of the series 
were therefore merely used as a check for the solution. We 
note that one using the Lyi only to fit this system would 
mistakingly derive a much higher H I (>20.0). Metals are 
found at the same redshift. 
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Figure 6. The two absorbers detected towards BRI J0137-4224 
at ^abs=3.101 (logAfHi=19.81±0.10) and 

^abs=3-665 (log AfHi=19-ll±0.10). The latter system is detected 
down to Lyll, but for display purpose only the first 7 members 
of the series are shown. 
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Figure 7. The absorbers detected towards PSS J0209+0517 at 
2,bs=3.666 (logArHi=20.47±0.10), 2abs=3.707 

(logArHi=20.55±0.10) and 2abs=3.863 (log AfHi=20.45±0.15). 
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Figure 8. The absorber detected towards BR J0311-1722 at 
2,bs=3.734 (logAfHi=19.48±0.10). 



(x) BR J0334-1612 (zcm=4.363): A DLA at 2abs=3.56 
was reported by Peroux et al. (2001) from medium reso- 
lution spectroscopy. They derive logA^Hi=21.0 in excellent 
agreement with the measurement made here from higher res- 
olution spectroscopy: log Ai'Hi=21.12±0.15 with 2:abs=3.557. 
Nevertheless, we note the asymmetry in the shape of the 
profile, a possible signature of multi-component absorbers. 
Metal lines associated with this system are detected in the 
red part of the spectrum. Figure 9 shows our best fit of this 
H I line. We also note in the spectrum of this quasar a real 
and significant broad dip in the continuum at A ~ 5200 A. 
This cannot be explained by any known broad absorption 
lines, nor by an ion at the emission redshift and does not 
correspond to the beginning of the Lyman/3 forest either. 
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Figure 9. The DLA absorber detected towards BR J0334-1612 
with 2abs=3.557 (log AfHi=21.12it0.15). We note the asymme- 
try in the shape of the profile, a possible signature of multi- 
component absorbers. 




-2000 



-1000 1000 

Relative Velocity (km/s) 



200C 



Figure 10. The sub-DLA absorber detected towards BR 
J0419-5716 at Zabs=3.063 (log ArHi=19.17±0.10). 



Similar features (although not as dramatic) have been re- 
ported by Levshakov et al. (2004). In that case, the dips are 
clearly associated with O vi at the emission redshift and are 
interpreted as ejected material or intracluster cooling flow. 

(xi) BR J0419-5716 (2cm=4.461): This quasar was ob- 
served by Peroux et al. (2001) who did not detect any 
sub-DLA in its medium resolution spectrum. Using higher- 
resolution data, we detect a sub-DLA with 2;abs=3.063 with 
log A'hi=19.17±0.10 as well as one metal transition. Fig- 
ure 10 shows our best flt of this H I line. 

(xii) BR J2017-4019 (2em=4.131): This quasar is a Broad 
Absorption Line (BAL) quasar (Peroux et al. 2001). We 
carefully search for DLA/sub-DLAs in the regions free from 
BAL features but do not detect any. We note however O vi 
AA1032 and 1037 features at Zabs=3.9966, 3.9692 and maybe 
also at 2:abs=3.7300. 

(xiii) PSS J2154+0335 (2om=4.363): This quasar was ob- 
served by Peroux et al. (2001) who report a DLA at 
2abs=3.61. This system falls in the setting gap of our UVES 
spectra. On the other hand, we discover a new sub-DLA at 
2abs=3.177 with log A''hi = 19.23±0.15. No metals were found 
associated with system over the limited wavelength range of 
our spectrum. Figure 11 shows our best fit of this H I line. 
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Figure 11. The sub-DLA absorber detected towards PSS 
J2154-I-0335 at Zabs=3.177 (log AfHi=19.23±0.15). 
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et al. (1996b) who did not 
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Figure 12. The sub-DLA absorbers detected towards BR 
J2215-1611 at 2abs=3.656 (log AfHi=19.01±0.15) and 2abs=3.662 
(log AfHi=20.05±0.15) respectively. 



find any DLA in the spectrum. Here we report on two 
new sub-DLAs with log iVHi = 19.01±0.15 and 20.05±0.15 at 
-Zabs=3.656 and 3.662 respectively. The two are separated by 
just ~320 km s~^. The Lyman series down to Ly3 is detected 
where the division into two distinct systems is unambiguous. 
Several metal lines are detected at these redshifts too. Fig- 
ure 12 shows our best fit of these H I lines. We also note OVI 
AA1032 and 1037 features at Zabs=3.9785 along this quasar 
line of sight. 

(xv) BR J2216-6714 (zcm=4.469): This quasar was 
observed by Peroux et al. (2001) who do not report 
any absorber but mention a possible sub-DLA candidate 
at 2abs=3.37. We now confirm that system and derive 
log Ai'Hi=19.80±0.10 at 2abs=3.368. Several metal lines are 
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Figure 13. The sub-DLA absorber detected towards BR 
J2216-6714 at Zabs=3.368 (log ArHi=19.80±0.10). 
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Figure 14. The DLA absorber detected towards BR J2239-0552 
at 2abs=4.079 {logiVHi=20.55±0.10). 



also observed at this redshift. Figure 13 shows our best fit 
of this H I line. 

(xvi) BR J2239-0552 (2cm=4.558): This quasar was ob- 
served by Storrie-Lombardi et al. (1996b) who report a 
DLA at 2abs=4.08. We now confirm that system and de- 
rive log Ai'Hi=20.55±0.10 at Zabs=4.079 in good agreement 
with estimates from medium resolution spectroscopy. The 
absorber is detected down to the LyS and has associated 
metal lines. Figure 14 shows our best fit of these H I lines. 

(xvii) BR J2349-3712 (2em=4.208): This quasar was ob- 
served by Peroux et al. (2001) at medium resolution. They 
do not find any DLA along its line of sight. Here, we re- 
port on the discovery of three sub-DLAs, two of which 
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Figure 15. The sub-DLA absorbers detected towards BR 
J2349-3712 at Zabs=3.581 (log AfHi=19.12±0.10) while the two 
others are at ^abs— 3. 690 (logArHi=19.79±0.15) and Zabs=3.696 
(logArHi=19.78±0.10). 



are blended. The first system has log A^hi=19.12±0.10 at 
-2abs=3.581 while the two others have log AfHi=19.79±0.15 
at 2abs=3.690 and log ArHi=19.78±0.10 at 2abs=3.696 (i.e. 
A.V ~380km s~^). The subdivision of these is unambiguous 
in Ly2. For all the systems the Lyman series is available 
down to Ly3. No associated metal lines are detected over 
the limited wavelength coverage of our spectrum. Figure 15 
shows our best fit of these H I lines. 
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Table 3. This tabic compiles the DLA column density and red- 
shift estimates from 5A (FWHM), signal-to- noise ratio per pixel 
~ 20 (Pcroux ct al. 2001; except BR J2239-0552 which is from 
Storrie-Lombardi et al. 1996b) and 2A, signal-to-noise ratio per 
pixel ~ 25 quasar spectra from the present study. 



Quasar 



BR J0006-6208 
PSS J0133H-0400 

PSS J0209-f0517 

BR J0334-1612 
BR .12239-0552 



mean 
min value 
max value 



2.4 Medium/High- Resolution A'hi Comparison 

Seven of the quasar absorbers described above are classi- 
cal DLAs with logA/Hi>20.3. In Table 3, the Nm column 
density estimates from medium resolution spectroscopy of 
these DLAs is compared with the new high-resolution mea- 
surements presented here. The last two systems (toward BR 
J0334-1612 and BR J2239-0552) were also observed by 
Storrie-Lombardi & Wolfe (2000) who found A^hi in very 
good agreements with ours. The comparison shows that pre- 
vious A^hi estimates from SA resolution quasar spectra are 
reliable. In fact, it is also known that estimates of systems 
with high column density from echelle data might be affected 
by the difficulty to trace the correct continuum over differ- 
ent echelle order. In the sub-DLA regime however, high- 
resolution spectra are definitely required in order to accu- 
rately measure AThi. 



3 ANALYSIS 
3.1 Survey's Properties 

Survey sensitivity: Table 4 lists the minimum (zmin) and 
maximum («max) redshifts along which a sub-DLA could be 
detected along each quasar line-of-sight. Zmin corresponds 
to the point below which the signal-to-noisc ratio was too 
low to find absorption features at the sub-DLA threshold of 
Wrest = 2.5 A, and Zmax is 3000 km s~^ blueward of the Lya 
emission of the quasar. We took care to exclude the DLA 
regions and the gaps in the spectrum due to non-overlapping 
UVES settings when computing the rcdshift path surveyed. 

Figure 16 shows the cumulative number of lines of sight 
along which a sub-DLA could have been detected at the 5cr 
confidence level. This survey sensitivity, g{z), is defined by: 

g{z) = ^ i/(^r" - z)H{z - zrn (i) 

where H is the Heaviside step function. In Figure 16, it is 
compared with those of previous sub-DLA survey (Peroux 
et al. 2003b). It shows that our new observations probe a 
higher redshift interval than the first archive-based work 



Table 4. Rcdshift path surveyed. Zj^in corresponds to the point 
below which no flux is observed and 3000 km s-i blue- 

wards of the Ly-a emission line. Gap in non-overlapping settings 
and known DLAs are taken into account. 



Quasar 


Zcm 


^min 


^max 


BR J0006-6208 


4.455 


3.083 


3.188 






3.216 


3.732 






3.798 


4.450 


PSS JOl 18+0320 


4.230 


2.767 


4.225 


PSS J012H-0347 


4.127 


2.739 


3.410 






3.474 


4.122 


SDSS ,10124+0044 


3.840 


2.414 


3.835 


PSS J0133+0400 


4.154 


2.414 


3.247 






3.305 


3.676 






3.707 


3.765 






3.786 


4.149 


BRI J0137-4224 


3.970 


2.344 


3.965 


PSS J0209+0517 


4.174 


2.767 


3.658 






3.675 


3.855 






3.872 


4.169 


BRI J0244-0134 


4.053 


2.414 


3.247 






3.305 


4.048 


BR J0311-1722 


4.039 


2.560 


3.410 






3.474 


4.034 


BR J0334-1612 


4.363 


2.748 


3.541 






3.573 


4.358 


BR J0419-5716 


4.461 


2.748 


3.572 






3.639 


4.456 


BR J2017-4019 


4.131 


2.607 


2.994 






3.022 


3.407 






3.473 


4.126 


PSS J2154+0335 


4.363 


2.748 


3.572 






3.639 


4.358 


BR J2215-1611 


3.990 


2.315 


3.408 






3.473 


3.985 


BR J2216-6714 


4.469 


2.748 


4.464 


BR J2239-0552 


4.558 


2.894 


3.735 






3.801 


4.070 






4.096 


4.552 


BR .72319-3712 


4.208 


2.601 


3.107 






:-!. it:-! 


1.203 



and that the combination of the two provide a rather ho- 
mogeneous survey coverage in the range z=1.5 to 4.5. 

Redshift distribution: The sub-DLA sample described in the 
section above leads to 21 sub-DLAs along 17 quasar lincs-of- 
sight whilst we had 10 sub-DLAs making up our "statistical 
sample" toward 22 lower redshift quasars in our previous 
archival study (Dessauges-Zavadsky et al. 2003). We also 
confirm the presence of 7 DLAs along the same lines of sight. 
The histograms showing the rcdshift distribution of these 
homogeneous sub-DLA samples are shown in Figure 17. 

3.2 Sub-DLAs Statistical Properties 

Redshift Number Density: The number of quasar absorbers 
per unit redshift, n{z), is a direct observable. This quantity, 
however, can be used to constrain the evolution or lack of it 
only when deconvolved from the effect of cosmology. 

The data acquired here used in combination with recent 



medium resolution high resolution A A''hi 







•^abs 




(med-high) 


3.20 


20.9 


3.202 


20.80 


+0.10 


3.69 


20.4 


3.692 


20.68 


+0.28 


3.77 


20.5 


3.773 


20.42 


-0.08 


3.66 


20.3 


3.666 


20.47 


+0.17 


3.86 


20.6 


3.863 


20.43 


-0.17 


3.56 


21.0 


3.557 


21.12 


+0.12 


4.08 


20.4 


4.079 


20.55 


+0.15 


3.69 


20.58 


3.689 


20.64 


+0.08 


3.20 


20.30 


3.202 


20.42 


-0.17 


4.08 


21.00 


4.079 


21.12 


+0.28 
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Table 5. This table summarises the observed number density of quasar absorbers for different column density ranges together with the 
predicted number density of sub-DLAs (from Peroux et al. 2003a and Rao, Turnshek &; Nestor 2005; except for the sub-DLA number 
density which is from the present work). The empty entries are regions of redshift/column density parameter space as yet unobserved. 



log A^Hi 






> 17.2 












> 19.0p,od 






> 20.3 




z range 


# 


<z> 


dz 


n(z) 


# 


<z> 


dz 


n(z) 


n(z) 


# 


<z> 


dz 


n{z) 


0.01-1.78 


7 


0.64 


10.9 


0.64±0.24 










0.21 


n/a 


0.90 


n/a 


0.10±0.02 


1.78-2.50 










3 


2.06 


8.8 


0.52±0.20 


0.40 


34 


2.16 


186.1 


0.18±0.03 


2.50-3.00 


17 


2.82 


10.2 


1.67±0.41 


5 


2.74 


9.1 


0.73±0.25 


0.46 


15 


2.73 


84.5 


0.18±0.05 


3.00-3.50 


11 


3.14 


4.3 


2.57±0.77 


7 


3.15 


9.4 


1.01±0.28 


0.67 


15 


3.27 


57.0 


0.26±0.07 


3.50-4.00 


9 


3.84 


2.3 


3.89±1.30 


12 


3.74 


9.6 


1.62±0.36 


1.52 


16 


3.77 


42.4 


0.38±0.09 


4.00-5.00 


26 


4.29 


7.2 


3.63±0.71 


4 


4.13 


5.1 


1.28±0.40 


2.00 


9 


4.19 


18.5 


0.49±0.16 


total 


71 




34.9 




31 




42.0 






100 




514.8 





#: number of absorption systems. 

19.0*1^^,: #, <z> and dz refer to systems with 19.0<log AfHi<20.3, whilst n(z) is for all systems with 
log7VHi>19.0. 

n/a: "not applicable" refers to low -redshift Mg Il-selected systems and so the number of absorption systems 
and corresponding dz are not directly comparable to higher-redshift statistics. 
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Figure 16. Survey sensitivity function. The g(z) function shows 
the cumulative number of lines of sight along which a sub-DLA 
system could be detected. The light grey histogram represents 
our archival search (Peroux et al. 2003b) while the present work 
is represented in black. It shows that our new observations probe 
a higher redshift interval than the archival work. The dip at z ~ 
3.45 corresponds to the 5 quasars for which the gap in the 540 
setting was not covered. 



results from the literature allow us to determine this quan- 
tity for various classes of quasar absorbers. This is shown in 
Figure 18 where the number density for DLAs (Peroux et 
al. 2003a; Rao, Turnshek & Nestor 2005), sub-DLAs (both 
predictions from Peroux et al. 2003a computation and new 
direct measurements from the present work) and LLS are 
presented (also from Peroux et al. 2003a). It can already 
be seen that the predictions overestimated the number of 
sub-DLAs at z > 4 while they underestimated the number 
of such systems at z < 3.5. The observations show that the 
number density of sub-DLAs is flatter than expected. All 
these observations are tabulated in Table 5. Assuming no 
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Absorption Redshift 



Figure 17. Histogram of the redshift distribution of the sub- 
DLAs discovered as part of this survey. Peroux et al. 2003b refers 
to the statistical sample of the archival search undertook at lower 
redshift (grey histogram) . The black histogram represents the new 
high-redshift sample from the present study. 



evolution in the number density, $ and gas cross-section a, 
the lack of evolution in the redshift number density in a 
non-zero A— Universe can be expressed as (see for example 
Peroux et al. 2004a): 



n{z) — no(l + z) X 



H{z) 



Ho 



where 

H{z) 
Ho 



[Qmz{1 + zf - nA[z{z + 2)] + (1 + zf] 



1/2 



(2) 



(3) 



These no evolution curves are shown for each class of 
quasar absorbers in Figure 18 where no is taken to be n{z — 
0) = 0.045 ± 0.006 derived by Zwaan et al. (2005b) from 
high-resolution 21-cm emission line observations of the z — 
analogs of DLAs. The other curves are scaled according to 
a factor representative of the difference in redshift number 
density of sub-DLAs and LLS in the redshift range 2.5 < 
z < 3.0. Departure from no evolution at z > 3 is clear for all 
classes of quasar absorbers. For comparison with previous 
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Figure 18. Number density of various classes of quasar absorbers 
as a function of redshift taken from Peroux et al. (2003a) and Rao, 
Turnshelc & Nestor (2005) except for sub-DLAs (tliis work). Tlie 
horizontal error bars are the bin sizes and the vertical error bars 
are the l-cr uncertainties. The dotted bins are the predictions of 
the number density of sub-DLAs from Peroux et al. (2003a), while 
the light grey bins at log7VHi>19.0 correspond to the observed 
number density presented in this paper. DLA's n{z = 0) data 
point is from the 21-cm emission line observations of Zwaan et 
al. (2005b). The curves represent a non-evolving population for 
a non-zero A— Universe except for the dashed line which is for a 
A = 0, (Jo = Universe. 



work, the dashed line shows the non-evolving population for 
a A = 0, go = Universe. 

Size of Quasar Absorbers: With the assumption that all 
classes of absorbers (i.e. DLAs, sub-DLAs and LLS) and all 
members of a given class are arising from the same under- 
lying parent population, one can estimate the cross-section 
radius of spherical absorbers, a, from the observed redshift 
number density n{z). Following Tytler (1981) and using a 
non-zero A-cosmology, we find: 

-1 



n{z) 



c 

Ho' 



1 + ^) 



H{2 



Ho 



e^{L)KTiR^{L)d{L) (4) 



where Q. m is the matter density, Q.a is the contribution of the 
cosmological constant, R{L) is the average H I absorption 
cross-section radius of a spherical galaxy with luminosity L, 
t is the fraction of galaxies which have gaseous absorbing 
envelopes with filling factor k. $(L) is estimated from the 
Schechter (1976) galaxy luminosity function: 

If we further assume that a Holmberg (1975) relation be- 
tween optical radius and luminosity holds at all redshifts 

R ( L 



R* \lJ 

we obtain the following relation for the radius: 
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Figure 19. The redshift evolution of the characteristic radius, 
ii» , of DLAs, sub-DLAs and LLS. For comparison, the impact pa- 
rameters of spectroscopically confirmed z < I DLAs (see Boissier, 
Peroux & Pettini 2003 for a compilation) and Mgll-selected galax- 
ies, a proxy for LLS, (Steidel et al. 2002) are also plotted. 



where r(a;) is the Gamma function. We take e=l, k=1 and 
t—OA (derived from Peterson et al. 1979). Two different sets 
of parameters are used for the luminosity function: $* — 
0.0149 ± 0.04 h^ Mpc"^ and s = 1.05 ± 0.01 for z < 0.75 
from the Sloan measurements of Blanton et al. (2003) and 
<!>* = 0.0142l° °°i^ h^ Mpc-=* and s = 0.501^6 ^r z > 0.75 
from infra-red selected galaxies of Chen et al. (2003). These 
yield a H l gas radius of: 



R* 



. . -1 ^(^) 

^"lOO 



1/2 



(1 + -) 



H{z) 



1 1/2 



Ho 



[kpc] 



(8) 



where A=76 for z < 0.75 and A=91 for z > 0.75. The red- 
shift evolution of the characteristic radius of DLAs, sub- 
DLAs and LLS are plotted in Figure 19. All classes of 
absorbers have their characteristic radius, decreasing 
with decreasing redshift, and this effect is stronger at lower 
column densities. For comparison, the impact parameters 
of spectroscopically confirmed z < 1 DLAs (see Boissier, 
Peroux & Pettini 2003 for a compilation) and Mgll-selected 
galaxies, a proxy for LLS, (Steidel et al. 2002) are also plot- 
ted. Given the number of assumptions made in the calcu- 
lation, the few available observations are in relatively good 
agreement with our results. 

Column Density Distribution: The differential column den- 
sity distribution describes the evolution of quasar absorbers 
as a function of column density and redshift. It is defined 



f{N,z)dNdX 



■dNdX 



(9) 



where n is the number of quasar absorbers observed in 
a column density bin [A'', A'^ + AA'^] obtained from the obser- 
vation of m quasar spectra with total absorption distance 
coverage ^^j^AXi. The column density distribution for 
two redshift ranges, z<3.5 and z>3.5 are shown in Figure 20. 
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Figure 20. Column density distributions for two redshift ranges down to the sub-DLA definition. The iiorizontal error bars are the bin 
sizes and the vertical error bars represent the uncertainties. The dotted bins are the predictions from fitting a Schechter function to the 
expected number of LLS (Peroux et al. 2003a), while the two solid bins at 19.0<log AfHi<20.3 correspond to the direct observations from 
the sample of sub-DLAs presented here. 



Table 6. This table summarises the observed column density distribution of quasar absorbers down to the sub-DLA column density. The 
empty entry in the highest column density bin at the highest redshift range reflects the lack of systems with log AfHi>21.35 at z>3.5. 



z range 






1.78- 


-3.5 








3.5- 


-5.0 




log A^Hi 




log 


log 


log 


log 




log 


log 


log 


log 


range 


# 


<Nm> 


f(N,z) 


fmin(N,z) 


fmaa;(N,z) 


# 


<Nm> 


f(N,z) 


fmin(N,z) 


fmaa;(N,z) 


19.00-19.50 


6 


19.23 


-20.5 


-20.4 


-20.8 


9 


19.22 


-20.2 


-20.0 


-20.3 


19.50-20.30 


9 


19.94 


-21.2 


-21.1 


-21.4 


7 


19.81 


-21.1 


-21.0 


-21.4 


20.30-20.65 


33 


20.45 


-21.9 


-21.9 


-22.0 


20 


20.44 


-21.6 


-21.5 


-21.7 


20.65-21.00 


20 


20.81 


-22.5 


-22.4 


-22.6 


6 


20.73 


-22.4 


-22.3 


-22.7 


21.00-21.35 


13 


21.08 


-23.0 


-22.9 


-23.2 


3 


21.07 


-23.1 


-22.9 


-23.5 


21.35-21.70 


5 


21.42 


-23.8 


-23.6 


-24.1 














These redshifts are chosen to allow a direct comparison with 
the work of Peroux et al. (2003a). In particular, the highest 
redshift is set to z—5 to match DLA surveys even though 
the sub-DLA search is made only up to z=4.5 (see Fig- 
ure 16). Again, the observations (solid bins) are compared 
with the predictions from Peroux et al. (2003a). The new 
data allow us to determine f{N,z) down to log A^hi = 19.0. 
The column density distribution has often been fitted with a 
simple power law (i.e. Prochaska & Herbert-Fort 2004), but 
there have been recent work showing that a Schechter-type 
of function is more appropriate (Pei & Fall 1995; Storrie- 
Lombardi, McMahon & Irwin 1996b; Peroux et al. 2003a). 
In any case, it should be noted that none of these two func- 
tions is physically motivated but rather are chosen so as to 
best describe the data. However, as more observations are 
available, a clear departure from the power law is observed. 
This fiattening of the distribution in the sub-DLA regime 
is indeed expected as the quasar absorbers become less self- 
shielded and part of their neutral gas is being ionised by 
incident UV flux. We note once more the paucity of very 
high column density DLAs at high-redshift. All the data are 
summarised in Table 6. 



3.3 JIhi+HcII Gas Mass 

Cumulative Number of Absorbers: In Figure 21, the number 
of quasar absorbers per unit comoving distance interval is 
plotted as a function of logA'ni- The new observations in- 
crease by one dex the A^hi column density parameter space 
probed and therefore provides a much better handle on the 
true shape of the distribution. This, in turns, is crucial to 
constrain the total H I gas content which is a function of 
both A^Hi column densities and redshifts. In our surveys, the 
DLA and sub-DLA samples (and associated redshift paths) 
are totally independent (Figure 21) in the sense that the 
quasars used to make up the samples are different in both 
cases given that higher resolution is required for the study 
of sub-DLAs. Therefore, the shape of both sections of the 
curves are unrelated and the fact that the function in the 
sub-DLA regime is a smooth continuation of the DLA part 
is resulting from the observations. 

f2Hi+Hcii.' The gas mass density, f2Hi+Hoii, observed in high- 
redshift quasar absorbers is classically expressed as a frac- 
tion of today's critical density; 
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Table 7. Redshift evolution of the gas mass density Ojji+HcII contained in classical DLAs and in sub-DLAs and expressed as a fraction 
of today's critical density. The total amount of H I gas is given in bold. 



log A^HI range 






log 
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dz 


dX 


xlO"^ sub-DLA 


quasar 


S^HI+HcII 


1.78-2.50 


2.17 


189.0 


586.0 


0.80±0.22 


36 


417 


2.06 


8.8 


27.6 


0.19±0.13 3 


35 


0.99±0.22 


2.50-3.00 


2.73 


87.2 


301.0 


1.02±0.31 


18 


274 


2.74 


9.1 
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0.17±0.10 5 


28 


1.18±0.31 


3.00-3.50 


3.25 


64.8 
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18 


178 


3.15 


9.4 


34.9 


0.22±0.11 7 


35 


1.32±0.31 


3.50-4.00 


3.77 


49.5 
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0.80±0.23 


19 
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3.74 


9.6 


37.9 


0.23±0.08 12 
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Figure 22. Redshift evolution of the gas mass density f^Hi+Holl expressed as a fraction of the critical density. The stellar mass density 
today is represented by f2* measured by Cole et al. (2001) for a Salpeter Initial Mass Function and the baryons in the hot gas is 
represented by Owhim measured from the incident of OVI absorbers by Danforth & ShuU (2005). The triangle at z=0 corresponds to 
the H I mass measured with radio observations in local galaxies (Zwaan et al. 2005a). The two error bars at z<2 are HHi+Hell from 
Mgll-selected DLAs (Rao, Turnshek Sz Nestor 2005). Left panel: the light grey bins at z>2 arc f2Hl+Hcll in systems with log A^hi>20.3, 
while the black bins correspond to the total H I gas, including the H I contained in sub-DLAs. Right panel: SIhi+HoII derived from the 
first Data Release of the Sloan survey is shown in light grey (Prochaska & Herbert-Fort 2004) together with our measurements of the 
total as mass (black bins). 



nHi{z) = 



AX 



(10) 



where fi = 1.3 is the mean molecular weight and mn 
is the hydrogen mass. The total gas mass, including H I gas 
in systems below the canonical DLA definition are plotted 
in Figure 22 and compared with the stellar mass density 
today (r2«) and H I mass measured with radio observations 
of local galaxies (triangle at z=0). Measurements at z<2, 
are nHi+Hcii from iVlgll-selected DLAs (Rao, Turnshek & 
Nestor 2005). The left panel of the figure decomposes the 
H I mass contained in both classical DLAs and in DLAs 
-I- sub-DLAs. These results are tabulated in Table 7. The 
right panel of the same figure shows f2Hi-i-Hoii derived from 
the homogeneous survey of DLAs found in the first Data 
Release (DRl) of the Sloan Sky Digital Survey (Prochaska 
& Herbert-Fort 2004). 



3.4 Clustering Properties 

As already noted in Peroux et al. 2003b, quite a few of the 
21 sub-DLAs which make up the new sample are close in 
redshift. The notes in the last column of table 2 emphasises 
in particular these systems less than ~400 km s~^ apart, 
which in practice means that they needed to be fitted to- 
gether: 2abs=3.995/3.999 toward PSS J0133-f0400 (Av = 
240 km s"^), 2abs=3.656/3.662 toward BR J2215-1611 
{Av = 380 km s'^) and Zabs=3. 690/3.696 toward BR 
J2349-3712 (also have An = 380 km s"^). Another example 
is the two Zabs=3.692/3.773 DLAs along PSS J0133-I-0400 
{Av = 5100 km s~^), with in between a 2abs=3.760 quasar 
absorber just below the sub-DLA definition log Ai'Hi< 19.0. 
In fact, the line of sight toward PSS J0133-f 0400 is partic- 
ularly rich, containing a total of 6 DLAs/sub-DLAs, all of 
which but one being separated by no more than Av = 15000 
km s~^. This complex group of systems appears to be alike 
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Figure 21. Cumulative number of quasar absorbers as a function 
of column density for two redshift intervals. Tiie discontinuity at 
log A'^HI=20.3 illustrate that the DLA and sub-DLA samples (and 
associated redshift paths) are made of different quasar spectra 
and are therefore totally independent. The observations show that 
the incidence of low column density absorbers is bigger at high 
redshift, as predicted by Peroux et al. (2003a). 

the multiple DLAs (MDLAs) studied in Lopez & Ellison 
(2003) and which are found to have abundance patterns 
distinctly different from the one of classical DLAs. In the 
present data, we do not cover the red part of the spectrum 
which would allow to detect metal lines to test further the 
hypothesis from Lopez & Ellison (2003). Acquiring a spec- 
trum covering the red part of this quasar would allow to test 
further the hypothesis about the possible truncated star for- 
mation of these MDLAs. 



4 DISCUSSION 

4.1 On the Ionized Fraction of Sub-DLAs 

The extension of the classical DLA definition to log A^'hi > 
19 proposed by Peroux et al. (2003a) has triggered concerns 
about the ionized fraction in these sub-DLAs. This issue is 
not relevant when quasar absorbers are counted up in order 
to quantify the total H I gas mass of the Universe, S7hi+hcII. 
Indeed, it is the sheer number of sub-DLAs which makes 
them add up to the classical DLA contribution, and this re- 
gardless of the amount of gas ionized in these systems. On 
the contrary, the quasar absorbers which have a high ion- 
ization fraction should not be included when measuring the 
neutral gas mass, as they trace hotter gas. Nevertheless, the 
evolution of neutral gas mass most probably is the prod- 
uct of several phenomena including star formation, galactic 
feedback, ionization of neutral gas as well as formation of 
molecular H2 gas (which will in turn leads to the formation 
of stars) . Which of these are the dominant processes probed 
by DLAs and/or sub-DLAs is at present unclear. 

When sub-DLAs are used for estimating the global 
metallicity, one might rightly wonder about the fraction of 
the gas (and consequently metals) in ionized form. In or- 



der to address this question, Dessauges-Zavadsky et al. 2003 
have carefully modelled every one of the 12 sub-DLAs mak- 
ing up their sample using detailed CLOUDY model. Overall, 
their findings demonstrate that indeed the ionized fraction 
X = /H account for 2/3"^^^ of the gas. But they have 
also shown that relative or absolute abundances estimated 
from X'^ /H° measurements are a good approximation of i.e. 
Xtot/ Fetot or Xtot/Htot- Therefore, using the metals com- 
monly detected in sub-DLAs for a measure of the metallicity 
of the H I gas in the Universe does not introduce a system- 
atic bias, even though the ionization fraction of the systems 
is not small. 

4.2 On the Nature of Sub-DLAs 

By assuming that both DLAs and sub-DLAs trace the same 
underlying parent population, a natural explanation for the 
nature of sub-DLAs could be that they are the outermost 
part of galaxies (Peroux et al. 2003b). This is illustrated 
by the absorber sizes calculations presented here where the 
characteristic radius of sub-DLAs is around 40 hj^pj, kpc and 
the one from DLAs is 20 hj~QQ kpc. 

The metallicity of sub-DLAs also seem to differ from 
the one of classical DLAs. Smoothed particle hydrodynam- 
ics simulations (Nagamine, Springel & Hernquist 2004a; 
Nagamine, Springel & Hernquist, 2004b) indicate that DLAs 
should be 1/3'^'* solar at z = 2.5 and even more metal-rich 
toward lower redshifts. Indeed there are lines of evidence 
pointing toward lower column density quasar absorbers like 
sub-DLAs being more metal rich at « < 2 (Peroux et al. 
2003b). This could be explained by classical DLAs being 
dustier than their sub-DLAs counterparts thus providing for 
the selection of against their background quasar in the first 
place. Recent computations have shown that even current 
radio-selected quasar samples (Ellison et al. 2001) are not in 
disagreement with such a scenario (Vladilo & Peroux 2005). 
In fact, from figure 9a of Vladilo & Peroux (2005), it can 
be appreciated that in the decade log A^hi=20— 21, from 10 
to 50% of the DLAs might be missed due to dust obscura- 
tion. Extrapolations of these mathematical formulations to 
the log Ai'Hi=19— 20 decade suggest that only less than 10% 
of the sub-DLAs might be missed by dust obscuration. 

If confirmed, this can be explained by the fact that in 
sub-DLAs, the Zn column density threshold does not com- 
bine with the HI threshold logA^Hi>20.3 that prevent their 
detection. 

4.3 Cosmological Evolution of the Mass Densities 

Figure 23 shows the cosmological evolution of some of the 
observable baryons in the Universe. In recent years, new 
observations have considerably changed the global picture 
with respect to previous studies (Lanzetta et al. 1991; Wolfe 
et al. 1995; Storrie-Lombardi et al. 1996a; Storrie-Lombardi 
& Wolfe 2000). In Figure 23, the direct comparison of the 
mass density in stars and in H I illustrate several issues: 

(i) S7hi+HoII observed in quasar absorbers at high redshift 
{z > 2) is low compare with the mass density observed in 
stars today, Q,*. 

(ii) fiHi+Heii observed in quasar absorbers evolve little 
over time from 2 = 5 to z > 0. 
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Figure 23. Observable baryons in the Universe as a function of 
time. The curve represents the mass density in stars from Rudnick 
et al. (2003) integrated from the Star Formation Rate measured 
by Cole et al. (2001). The uncertainties in that measurement are 
estimated to be around 15% but are also a strong function of the 
chosen Initial Mass Function. The error bars represent the mass 
density in the quasar absorbers deconvolved from the local critical 
density. A direct comparison of these two quantities illustrates the 
puzzling current situation concerning the baryon content of the 
Universe. 

(iii) SIhi+hcII observed in quasar absorbers at interme- 
diate and high redshift {z > 0) is high compare with H I 
inferred from 21-cm observations of local galaxies. 

(iv) flHi+Heii observed in quasar absorbers over most of 
the age of the Universe (since the Universe was 3 Gyrs old, 
corresponding to about z < 2) is not well constrained. 

The first point i) arises mostly from the fact that the 
most recent studies make use of modern cosmological param- 
eters (Spergel et al.2003), including a non-zero A cosmolog- 
ical constant, which compresses the high-redshift measure- 
ments with respect to the local ones. The possibility that 
large number of quasar absorbers are missing in optically 
selected quasar surveys is still hotly debated. While radio 
surveys looking for DLAs in quasar samples without optical 
limiting magnitudes (Ellison et al. 2001; Ellison et al. 2004, 
Prochaska, private communication) show that there is not a 
large number of DLAs missing, other observational indica- 
tors suggest that the DLAs we know of are not extremely 
dusty (Pettini et al. 1997; Ledoux, Petitjean & Srianand 
2003; Kulkarni et al. 2005; Murphy & Liske 2004). It should 
be emphasised however that there are two separate issues: 
i) what is the dust content of the quasar absorbers we know 
of today and ii) which fraction of the quasar absorbers are 
missed because their background quasar is not selected in 
the first place. An extension of Fall & Pei (1993) calculations 
taking into account the most recent observations addresses 
both these issues: Vladilo & Peroux (2005) show that whilst 
the dust content of the DLAs in current sample is not too 
high, the missing fraction is possibly quite important rang- 
ing from 30% to 50% at z = 1.8 - 3.0. 

Concerning points ii) and iii), the number of quasar 
absorbers observed at z > 2 is now reaching hundreds. The 



mild fluctuations in the redshift evolution in the range z = 5 
to 2 = 2 is within the small error estimates. Therefore the 
cosmological evolution of the total gas mass, IIhi+HcII, can 
be approximated to constant in that redshift range. Yet, 
radio observations of very large numbers of local H l-rich 
galaxies such as the HIPASS survey (Zwaan et al. 2005a) 
indicate a low value of fiHU-Hoii at z = 0. This would imply 
a fast evolution of the gas mass between z ~ 0.61 (Rao, 
Turnshek & Nestor 2005) and z = 0. 

The last point illustrates how challenging it still is to 
find and study quasar absorbers at intermediate redshifts. 
Warnings about f2Hi+Hoii measurements based on small 
number statistics have already been issued (Peroux et al. 
2004b) and unfortunately future developments in that direc- 
tion appear limited given the presently restricted availability 
of UV optimised instruments. 



5 CONCLUSION 

We have presented a new sample of high-redshift sub-DLAs 
{Nm> 10^^ cm~^) found in the spectra of 17 2 > 4 quasar 
spectra observed with the Ultraviolet- Visual Echelle Spec- 
trograph (UVES) on VLT. The statistical properties of this 
sample of 21 new sub-DLAs is analysed in combination with 
another 10 sub-DLAs from previous ESO archive studies. 
This homogeneous sample allows to determine redshift evo- 
lution of the number density of DLAs, sub-DLAs and com- 
pare it with that of LLS taken from the literature. All these 
systems seem to be evolving in the redshift range from z — 5 
to z ~ 3. Assuming that all the classes of absorbers arise 
from the same parent population, estimates of the char- 
acteristic radii are provided. Rt increases with decreasing 
column density, and decrease with cosmological time for all 
systems. The redshift evolution of the column density distri- 
bution, f(N,z), down to A^Hi= 10^^ cm~^ is also presented. 
A departure from the usually fitted power law is observed 
in the sub-DLA regime. f(N,z) is further used to determine 
the total H l gas mass in the Universe at z>2. The complete 
sample shows that sub-DLAs are important at all redshifts 
from z = 5 to z = 2 and that their contribution to the total 
gas mass Hhi+HcII is ~20% or more if compared with the 
Sloan results (Prochaska & Herbert-Fort 2004). Finally, we 
discuss the possibility that sub-DLAs are less affected by 
the effects of dust obscuration than classical DLAs. 
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